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Braneworld models with induced gravity have the potential to replace dark energy as the explana- 
tion for the current accelerating expansion of the Universe. The original model of Dvali, Gabadadze 
and Porrati (DGP) demonstrated the existence of a "self-accelerating" branch of background solu- 
tions, but suffered from the presence of ghosts. We present a new large class of braneworld models 
which generalize the DGP model. Our models have negative curvature in the bulk, allow a second 
brane, and have general brane tensions and localized curvature terms. We exhibit three different 
kinds of ghosts, associated to the graviton zero mode, the radion, and the longitudinal components 
\^ ' of massive graviton modes. The latter two species occur in the DGP model, for negative and positive 

I brane tension respectively. In our models, we find that the two kinds of DGP ghosts are tightly 

. correlated with each other, but are not always linked to the feature of self-acceleration. Our models 

' are a promising laboratory for understanding the origins and physical meaning of braneworld ghosts, 

' and perhaps for eliminating them altogether. 
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Models with extra spatial dimensions provide new perspectives on many long standing problems of particle physics 
and gravity. One of the most exciting suggestions is that an extra dimension, rather than dark energy, may be 



■ the origin of the currently accelerating expansion of the Universe. The original simple model of Dvali, Gabadadze 
and Porrati (DGP) posits that we live on a codimension one brane in an infinite flat five-dimensional bulk i. As 
generalized in ^^], the model has only two relevant input parameters: a crossover scale rc determined by the ratio of 

^-H ■ the 4d brane-localized gravity coupling with the 5d bulk one, and the tension of the brane. The Einstein equations 
^ : have two kinds of solutions: a normal branch and a "self-accelerating" branch; both branches cause a modification of 

■ the effective 4d Friedmann equation, with the second one having the property that at sufficiently late times it gives 
, an accelerated expansion without the need of dark energy to drive it. 

The DGP model has encountered many difficulties, but the one that looks most serious is the presence of ghosts 
CIhI in the weakly-coupled long distance regime of the theory linearized around the self-accelerating background solution 
^ ■ 11-0 ■ The physical origin and meaning of these ghosts are not yet clear. 

Ghosts, like tachyons, are often an indication of an instability of the theory perturbed around a given background. 
Indeed it is not surprising that gravitating systems with nontrivial backgrounds including branes (in the absence of 
supersymmetry) exhibit a variety of instabilities. In the case of DGP, it is natural to ask whether modified brane 
setups might avoid ghosts and tachyons, while still exhibiting the key feature of self-acceleration. 

In this paper we present a large class of warped brane setups that generalize the DGP model. Our basic idea is to use 
an AdS^ bulk space, taking advantage of the fact that slices of AdS^ can be flat, AdSi, or ^54. We consider models 
with one or two branes, taking the brane tensions and the brane-localized gravity couplings as input parameters, 
along with the bulk cosmological constant. We obtain background solutions with both normal and self-accelerating 
branches, for branes that are either ^^54 or ^54. The self-accelerating branches for the de Sitter branes give models 
precisely analogous to the DGP model; the self-accelerating branches for AdS4 branes are not of obvious cosmological 
interest, but allow us to generalize the well- understood physics of the Karch-Randall model 

In two previous papers [1, Q we set up the necessary tools to study generic warped models with arbitrary brane 
tensions and localized curvature terms. The size of the extra dimension and the brane curvature are well defined 
functions of the input brane parameters, and therefore by varying the latter we can parametrize a large number 
of different models. In one limit we make explicit contact with the original DGP model, while in another limit 
we reproduce the AdS5/AdS4 Karch-Randall model. In [1, we studied the phenomenology of AdS4 branes for 
background solutions on the standard branch, exhibiting a variety of ghosts and tachyons in parts of the parameter 
space of models. In this article we will extend that study to the self-accelerating branch and also consider the case 
in which the branes are dS^. We will pay particular attention to the presence of ghosts in the spectrum and how this 
is related to which branch we are considering. 

The outline of the paper is the following. In Section [ll] we describe some basic facts about ghosts in theories with 
gravity. In Section IIIII we review the main features of our new models, extending the results of 0, Q to the self- 
accelerating branch and dS^ branes. The spectrum, with particular emphasis on the presence of ghosts and tachyons. 
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is studied for the different cases in Section HVl Section |V] is devoted to our conclusions, and some technical details 
are relegated to appendices. 



II. GHOSTS AND GRAVITY 



With the ( — h++) metric signature, a canonical scalar field kinetic action in flat space is written 



/ 



2' "^'■^"''•^ 2 
Reversing the sign of the first term produces a ghost: 



d^x l^-rj'^-'d^q^d,<j) - ^m'cb'j . (2) 

As discussed e.g. in p^ . defining the ghost propagator with the usual Feynman ie prescription results in a nonunitary 
theory with negative norm states; defining the propagator with the opposite ie prescription preserves unitarity but 
induces a catastrophic instability from propagating negative energy states. 

Thus kinetic ghosts are similar to tachyons, in that they can be regarded as instabilities. Indeed some ghosts are 
also tachyons: this is the case in ([2]) for > 0, since the solutions of the wave equation will be exponentials rather 
than plane waves. As with tachyons, in some cases it may be possible to cure a ghost instability by perturbing around 
a different ground state. Such a shift can be thought of as the nonperturbative formation of a tachyon or ghost 
condensate; a tachyon condensate may produce a stable static ground state while a ghost condensate may produce a 
time dependent but ghost-free ground state [lli • 

In addition to ghosts and tachyons, one might worry that simple scalar field theories can suffer from other diseases. 
Suppose we add a derivative self- interaction term to 

d'x (^-irjf^-^d^^d^cj) - im202 + _^(^M^a^0a^0)2^ ^ (3) 

where A is the UV cutoff of this effective theory. For c < 0, this theory appears to have analyticity problems and 
exhibits superluminal modes in the massless limit, despite a lack of tachyonic instability [12]. However Jenkins and 
O'Connell have recently shown [l^ that "positivity constraints" of this type are equivalent to a no-tachyon condition 
in a (partial) UV completion. In other words, effective theories of the form ((3|) with c < arise from integrating out 
tachyons before curing the tachyonic instability. Thus kinetic ghost and tachyon instabilities appear to exhaust the 
physically distinct diseases of scalar field theories. 

This simple picture, however, is intrinsic to flat space propagators and flat space kinetic terms. Adding gravity to 
the picture complicates the discussion of ghosts, and has led to much debate in the literature. Once we add gravity, 
we must be on the lookout for new kinds of ghosts. In addition to scalar ghosts, the longitudinal mode of a massive 
graviton can be a ghost. As we will see, even a massless graviton can be a ghost. Because gravity is a nonlinear 
theory with a large local symmetry, we must also work harder to understand the physical significance of ghosts. 



A. de Sitter ghosts 

It was shown many years ago by Higuchi that the longitudinal modes of massive gravitons can be ghosts if we 
are in a de Sitter background. Letting 127i^ denote the constant de Sitter curvature, massive gravitons with mass in 
the range < < are ghosts. These correspond to nonunitary representations of the de Sitter group 5*0(4, 1), 
as discussed in [l^, [l^, but the analysis of Deser and Waldron fl7] shows that they are also kinetic ghosts. As we 
will see in the next section using an explicit on-shell tensor decomposition, the effective kinetic term of a longitudinal 
massive graviton mode s(x) in a dS/^ background is given by: 

where is the de Sitter Laplacian = 9f//(t)^, with /(t) = exp(7^i) coming from the ^5*4 background metric: 
(7oo = —1, = f^Sij. From ^ we see immediately that the longitudinal mode is a ghost for < < 2T-p. It would 
also naively appear that the longitudinal mode is a tachyon for < < |7i^. However the analysis of tachyonic 
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instabilities in dS^ is complicated by the fact that the naive Hamiltonian corresponding to the wave operator in (j4|) 
has an explicit time-dependence from f(t) and is therefore not conserved. The operator corresponding to the actual 
conserved energy is manifestly positive inside the de Sitter horizon for any m? > 0, as shown in p^ . This is referred 
to as a "mild" tachyonic instability in the last reference in for our purposes we will not call such modes tachyons, 
reserving that name for cases such as m? < in ([4]) which resemble tachyons in flat space. 

It was shown in ^4] that the DGP model with positive brane tension contains a longitudinal massive graviton mode 
ghost on the self-accelerating branch. Intriguingly, for negative brane tension the massive graviton ghost is absent, 
replaced instead by a ghost radion scalar. In the warped models presented in the paper, we will find that these two 
kinds of ghosts are also tightly related. 

The physical interpretation of ghosts in theories with gravity can be obscured by the presence of local symmetries. 
For example the DGP model in the limit of a tensionless brane has a massive graviton mode with m? = 2H^; 
the resulting 4d massive gravity theory has an enhanced "accidental" local symmetry beyond 4d general covariance 
[isl [Toj . Charmousis et al. found a ghost in this limit as well [3] , but Deffayet et al have argued ^] that this mode is 
actually a Lagrange multiplier enforcing an explicit gauge-fixing of the extra symmetry. This is an interesting special 
case to be considered elsewhere. 



B. when a ghost is not a ghost 

It can be argued (incorrectly) that gravity renders the physical meaning of ghosts ambiguous. Consider, for example, 
pure 4d Einstein gravity. Perform a conformal rescaling of the metric by the field redefinition 



9^i•^ ^ ( 1 + ) 5/.^ ' (5) 



where (jii^x) is a 4d scalar field and M is the reduced Planck mass. Now substitute the rescaled metric into the 
Einstein-Hilbert action, assuming that the background metric before rescaling was flat: 

3 1 

2 1 + 0/M 

This looks like a kinetic ghost. Now suppose we start with a Lagrangian whose perturbation theory around flat space 
is obviously ghost-free: 



-g M^R ^ - ^^7^ il''''d^cj)d,c^ . (6) 



(7) 



It appears naively that by a simple field redefinition of the metric we can produce a ghost. 

The flaw in this argument is that we have ignored the subtleties of general covariance and gauge- fixing (20i] . Consider 
again 4d Einstein gravity expanded around flat space: g^,y = rj^^i, + /i^j,. As discussed in ^g], the massless symmetric 
tensor /i^^ can be decomposed as 

= + df_,v^ + d^v^ + d^d^ipi + c^^ + d^n^ + d^n^ + 7]^„ip2 ■ (8) 

Here is traceless, transverse, and orthogonal to ti^, giving 2 degrees of freedom. Similarly is transverse and 
orthogonal to n^, giving 2 degrees of freedom, while the longitudinal null vector gives 1 more degree of freedom. 
Lastly, c^i/ is traceless but not transverse, giving 3 degrees of freedom, which together with the two scalars ipi and ip2 
adds up to the total 10 degrees of freedom of h^^^. 

Obviously v^, and ipi are the 4 pure gauge modes of 4d general covariance, while (3^i, gives the two traceless 
transverse propagating degrees of freedom of an on-shell massless graviton. The four remaining degrees of freedom 
represented by ip2 and c^i, are subject to constraints from the equations of motion. In the absence of sources they are 
constrained to vanish. Just as for the more familiar case of the time-like component of an abelian gauge field, these 
modes do not really propagate. In the post-Newtonian approximation [2T|, tp2 is proportional to the Newtonian scalar 
potential (analogous to the electrostatic potential) while c^^ contains the Newtonian vector potential. So ip2 represents 
the same kind of ghost as the time-like component of an abelian gauge field (which in the nonrelativistic limit gives 
the electrostatic potential). Such ghosts are not really ghosts since they do not propagate: the corresponding degrees 
of freedom are completely fixed in terms of sources by solving constraints from the equations of motion. 

Thus the metric rescaling JS]) in the theory described by (O mixes a non-ghost matter scalar with a pseudo-ghost 
mode of the metric. The full quadratic action of the rescaled theory will have kinetic terms mixing (j) with ip2. 
Presented with such an action ab initio, the proper way to extract its physical degrees of freedom is to gauge-fix 
it and diagonalize the kinetic terms. The only way to diagonalize the kinetic terms is by shifting ^2 ^ ^2 ^ 'i'l^l 
effectively undoing the metric rescaling ([5]). The theory is then obviously free of ghosts. 
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C. summary 

The moral of this story is that there is an unambiguous procedure for determining the presence of ghosts in a theory 
with gravity expanded around a given background solution. First, determine the local symmetries and gauge-fix them. 
Second, diagonalize the full quadratic action. Third, extract the propagators of the physical degrees of freedom and 
check for kinetic ghosts. It is also important to check for tachyons; in our ^5*4 example above, the longitudinal graviton 
mode s{x) for < is an example of a tachyon which is not a ghost. 

As we argued above, for scalars in flat space ghosts and tachyons exhaust the list of physically distinct pathologies 
for weakly coupled theories. This is less obvious for gravity. However if we avoid regions of strong coupling and 
abstain from integrating out degrees of freedom, there does not seem to be any reason to impose additional positivity 
constraints. We will assume this for our analysis. 

There is one other possible source of ambiguity. Because gravity is a highly nonlinear theory, one could imagine that 
there are physical setups in which even the long-range solutions are intrinsically nonlinear. In such a case ghosts of 
the linearized solutions are irrelevant, since the linearized solution is not itself an approximation to the real solution, 
even very far away from any source. It has been argued that this may occur in the DGP model Note that 

the claim that nonlinearities of gravity are important at cosmological scales is different from the claim that matter 
nonlinearities are important at cosmological scales |23| : the latter has also been suggested as a replacement for dark 
energy in explaining late time accelerated expansion [24|. At any rate we will not pursue this idea here, since we find 
the linearized analysis of our models quite challenging enough! 

III. WARPED MODELS WITH INDUCED GRAVITY 

In this section we will review the main features of the two-brane models introduced in extending the results 

to include the self-accelerating branch and the case of positively curved branes. The model is described by the 
following action 




+ [ d*xx/-g(''){2Mf'R.'-'^ - V,) + UI^ I K . (9) 
i J y=yi J OM 

This action represents a general warped gravity setup with codimension one branes, written in the interval picture. 
We have only displayed one interval since physics on the second interval is fixed by a Z2 symmetry. In ([9]) M is the 5d 
Planck scale, A — — 24M^fc^ is the bulk cosmological constant giving a bulk curvature k, the Mi are the coefficients 
of brane-localized curvatures "^^^^ the Vi are brane tensions and K is the trace of extrinsic curvature. As described 
in we will use the straight gauge formalism to keep the brane locations fixed at y — 0^ and y = L~ , even in the 
presence of linearized fluctuations of the metric. There are no "brane-bending" modes in a straight gauge. 
The background solution can be written as 

ds^ = GMNdx^'dx^ = g^.^dx^'dx" + (^?/^ (10) 

where y is restricted to the interval < y < L and 

g^iv ^ a^{y)it,v , (11) 

with 7^1/ the metric of dS^ or AdSi with constant curvature 127i^. The explicit form of the warp factor, a{y)^ the 4d 
curvature constant H^, and the coordinate size of the extra dimension L, depend on the input brane parameters that 
also determine how the AdS^ bulk is sliced. It is convenient to define the brane parameters in terms of dimensionless 
quantities, 

V, = kMf/M^, = V^/2kM^ . (12) 
The derived features of the model are determined by the following combinations of input parameters: 

T;^^^ (^-l±/l + W6+^) . (13) 

The choice of sign in (I13|) corresponds to the two independent branches (per brane) of the background solutions to 
the equations of motion. 
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A. slicing the AdSs bulk 



For a given set of five input parameters fc, wq, wq, vl, wl, the first question is to determine the possible ^5*4 or 
^^5*4 slicings of the AdS^ bulk. This is set by whether the absolute value of is larger or smaller than one. We show 
in Figured] the regions in the {v, w) plane that correspond to |T| > 1 with dark shade and |T| < 1 in light shade. In 
the unshaded regions T is complex. The left and right panels are for T"*" and T~ , respectively. The curved boundaries 
correspond to the line w ~ — 6(l + w^)/w, the horizontal lines that separate the dark and light areas correspond to 
w = ±12 (and v smaller or larger than ±1, depending on the panel) and finally the slanted line on the left panel 
separating the two light shades (with T+ positive or negative) corresponds to w = ~6v. 

Simply put, in the dark shaded regions the branes are both dS^, in the light shaded regions they are both AdS^. 
We parametrize the warp factor differently in each case: 




V 



FIG. 1: Regions of the {v,w) plane corresponding to different background solutions and signs of T. The left panel is for the 
r+ branch and the right panel is for the T~ branch. The two (lighter) darker shades correspond to {A)dS4 branes. In the 
unshaded regions the branes are neither AdS4 nor dS'4. 



• AdS5/AdS4: If \Ti\ < 1, then the branes are AdSi. The warp factor reads 

shk{y - : 
cosh kyo 



a{y) = (14) 



with the turning point yo and L given by 

To = tanhfcyo, Tl = tanhfc(L — yo). (15) 
Note that the requirement L > implies 

To+Tl> 0. (16) 
Finally the brane curvature is given by {H^ = —H^ < 0) 

cosh kyo ^ ^ 

• AdS5/dS4: If \Ti\ > 1, then the branes are dS'4. The warp factor reads, 

sinhk{y~yo) , . 

a{y) = r-r-T , (18) 

smh kyo 

with yo and L given by 

To = cothfcyo, Tl = cothfc(L - yo)- (19) 

Note that if < yo < the warp factor vanishes at some point in the bulk. This is a true singularity and 
therefore does not give rise to an acceptable model. Thus we have to require either yo < or yo > L for the 
case of ^5*4 branes. This requirement, together with L > 0, implies the conditions 

To • Ti < 0, To + Tl< 0. (20) 
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Finally the brane curvature is given by 



where eo is the sign of Tq . 



k 

^ ~^ sinh kyo ' ^ 



B. standard branches and self— accelerating branches 

For fixed input brane parameters, we have two different background solutions per brane, denoted by the sign choice 
Ci = ± in T/' (not to be confused with ii, the sign of itself). These two branches exhibit different physical features. 
In the following, we will focus our attention to the brane at y = 0. The restriction L > 0, in addition to j/o < or 
L < yo for dSi branes, will then imply restrictions on the parameters of the brane located at y — L, such that not all 
combinations of signs give sensible solutions. 

The choice of branch for the brane at y = 0, i.e. the choice of eo, corresponds to the standard and self-accelerating 
branches for the cosmology of brane-induced gravity. In our notation, the standard branch corresponds to Tq' whereas 
the Tq" branch is the self-accelerating one. The reason is the following. For fixed values of the brane parameters, 
Wo, Wo, the brane curvature is larger for the self-accelerating branch than for the standard one^. This is obvious for 
\wo\ < 12, for which corresponds to ^^5*4 solutions whereas Tq~ gives a ^5*4 brane (see Fig.[T]). In the rest of the 
parameter space, where we have either two ^^5*4 or two ^5*4 solutions, it is still true that 

n^[T+ivo,wo)] <n^[TQ{vo,wo)]. (22) 

Furthermore, for positive vq, the Randall-Sundrum [2^ tuning wq — +12 results in a flat brane (independently of the 
particular value of t^o) for and an inflating brane for Tg~, thus the names standard and self-accelerating branches. 
This is the exact warped analogue of the DGP tensionless brane. In the following we will loosely use the term branch 
to denote the choice of sign in either brane. However, as we said, it is the brane at y = that determines the standard 
or self-accelerating nature of the solution, and therefore we will only apply these terms to the brane at y ^ 0. 



C. properties of the spectrum 

The properties of the physical spectrum were described for the AdS4 brane case in [1, Q , and can be easily extended 
to the dS^ case. The spectrum is obtained by computing the action at the quadratic level for small perturbations 
around the background metric, 

Gmn Gmn + hMN- (23) 

This can be done in a straight gauge, for which the branes are straight and at flxed positions, y = 0~^,L~, and 
hp,i{x,y) = 0. A further gauge-fixing exhibits the radion, 

/i44(x,y) -F(y)7/.(x), (24) 

where F{y) is a fixed but arbitrary function (part of the gauge choice) and il>{x) is a four-dimcnsional scalar. The 
components of h^^{x,y) can then be shown to have the following 4d physical degrees of freedom: 

• massless modes: there is a graviton zero mode, B^^{x) (2 degrees of freedom) and the above mentioned radion, 

il}{x). 

h^,. ^ a^iy)B^,ix) + a^yiiy)V + g^My)^ . (25) 
where the y dependence in the radion piece is given by 

yi{y) = x:^iz ^ If ^ ^ , y2iy) ^ xk{z - If + -:f . (26) 



The name self-accelerating is actually a bit of abuse of language in these more general configurations as there are regions of parameter 
space where the Tq branch gives an AdS4 background, which is never accelerating. Nevertheless even in this case there is a net positive 
contribution to the brane curvature as compared with the standard branch. 
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In this last equation we have defined the variable z 



tanhfc(y - yo), < 0, 
cothfc(y ~ yo), > 0. 



(27) 



Note that we have Ti = 9iz{yi), with —Oq = 6l = 1. x-, ^ ^-nd 5^ are remnants of the residual gauge freedom 
(related to F(y)) whose detailed form is not important for the discussion below. 



• massive modes: there is a tower of massive gravitons (5 degrees of freedom each) 



(28) 



where q labels the mass, q = m? /Ti"^ , and 



,{q) ^(-l + v'5^)/2V 
^0 



(-l+V^^^)/2 



(g) ^(-l + v'Tr^)/2 



(29) 



where the P's and Q's are associated Legendre functions. The mass spectrum of modes is determined by solving 
the determinant equation, 



with 



(30) 



a* 
b, 



{VMT^ - 1) + mP:l+^^^,iO^T^) , 



20.(1 -I;^)^g^ 



z=eiTi 



(31) 
(32) 



All these equations are valid for both signs of the brane curvature, provided the right definition of z, the warp factor, 
a{y), and the sign of the brane curvature TP are used.^ Performing the integration over the extra dimension in the 
quadratic terms for the different modes we obtain the corresponding kinetic coefhcients (normalization constants). 
For the massless modes: 



^3 ^ 



k 
2k 



- kyi 



E 



Vj + Tj 
1 + v,T, 



(33) 
(34) 



and for the massive modes: 



2 / y'^^'dz + J2[v,{l~z')y(^^%^eT, 



-To 



(35) 



If any of these coefficients becomes infinite, the corresponding mode is not normalizable and decouples from the 
spectrum. If it vanishes, then we are in a region of strong coupling. If it becomes negative, then the mode is a ghost. 



^ The price for this convenient notation is that q as defined here is —q as defined in Q. 
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D. longitudinal graviton modes in dS4 

As already mentioned, for ^5*4 the kinetic terms of some massive graviton modes have an intrinsic mass-dependent 
kinetic coefficient aheady at the 4d level, in addition to the overall coefficient Cg"*' that they inherit from the Kaluza- 
Klein (KK) decomposition. Thus in this case we need to compute this intrinsic coefficient before deciding if the mode 
is a ghost. 

To compute these coefRcients, we need an explicit orthogonal decomposition of the five "helicities" of each massive 
4d graviton bI^J{x). The appropriate orthogonal decomposition turns out to be: 

B^^\x) = f'/'vU'^\x) + /-2ao9.(/i/2,(.)(^)) , (36) 
Blfix) ^ f'^'b<^{x) + r'do{f^'id,vfix) + d,v'i'^\x))) 

where 

Pv - - ^ (37) 

is a transverse projection operator in the fiat 3d sense, and Vf = df / f{tY is the 3d Laplacian defined already in the 
previous section. In ([55]) the h'f^ {x) are traceless-transverse in the flat 3d sense (2 degrees of freedom), corresponding 

to the helicity ±2 modes of the massive gravitons in the KK tower. The v'f\x) are transverse vectors in the fiat 3d 
sense (2 degrees of freedom), corresponding to the helicity ±1 modes of the massive gravitons in the KK tower. The 
s'^'^\x) are the longitudinal modes of the massive gravitons. 

Substituting the ansatz ([5^ . it is straightforward but tedious to show that the full 5d bulk equations of motion 
reduce to a set of constraint equations which are automatically satisfied, together with the following dynamical 4d 
equations of motion in the de Sitter background: 

72 _ (^-2 _ 5a^2^\ Aq), 



= (-a„^+v^-K-^7^^)j5,^J^(x), 

2 i _«2 , Y72 _ (2 _ 9a-2N\ Jq) 



4 



V -ao^ + v^-K-^7^^)urw, (38) 



Solutions to the equation of motion for the vectors v^^\x) and scalars s^'^^x) appear to be defined only up to an 
arbitrary harmonic function annihilated by Vf , but a residual 4d general coordinate invariance removes this ambiguity. 

We can also substitute our orthogonal decomposition ((36)) into the effective 4d quadratic action, to read off the 
intrinsic mass-dependent kinetic coefficients. Since the decomposition involves derivatives, this is really only mean- 
ingful with the derivatives evaluated on-shell. The interesting case is for the longitudinal modes, where after another 
tedious calculation we obtain the effective action shown in Q for each longitudinal graviton mode in the KK tower. 

Thus, even for values oi q — m? /Ti,^ such that Cg'"* > 0, the longitudinal graviton mode will be a ghost if < g < 2. 



IV. GHOSTS IN MODELS WITH INDUCED GRAVITY 



In this section we will discuss the spectrum of our models in the search of regions that are free of ghost and tachyonic 
instabilities. This is done by studying the solutions of the eigenvalue equation for the massive modes (|30p as well as 
the kinetic coefficients of these and the graviton zero mode and radion, p3ll35p . It is useful to separate the discussion 
for AdSi^ and dS^ branes, because of the different behavior of the longitudinal component of massive gravitons in the 
latter. Thus, in an ^^54 space, we only need to look at the kinetic coefficients of the different modes, ghosts being 
uniquely determined by the negative sign of their kinetic terms, whereas in ^5*4 space, the longitudinal component 
of a massive graviton with positive overall kinetic coefficient but with mass < < 2H^ is a ghost, becoming a 
(non-ghost) tachyon for < 0. 
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A. negative curvature: AdS^ branes 



The presence of ghosts for AdS^ branes was discussed, for the standard branch, in 9]. In this section we wiU review 
the results and generahze them to the self-accelerating branch. Let us start discussing the simpler case of the graviton 
zero mode. Recalling the definition of in the AdS4 background ([T5)) . we can write the kinetic coefficient for the 
graviton zero mode as 



1 t2 



tanh^^Tj 



T, 



1 - T/ 



(39) 



For AdS4 the global term outside the sum is positive so we only need to look at the sign of the sum. Let us define 
each term in the sum as 

Cf-^tanh-ir. + ^. (40) 

Cg°^'* is a growing function of Wi, for fixed Wj. Let us now consider its value at the boundaries of the ^^5*4 region 
in the {v,w) plane (light shaded areas in Fig. [1]). In the limit Vi — > ±cxd along the curve Wi = —6(1 + vf)/vi we get 

4°^'^ = ±oo. As we move along that curve towards the points Vi ±1 (and Wi = Tl2, respectively), it goes to =Foo. 
This is independent of which branch we have chosen for the branes. The rest of the boundaries, however, depend on 
whether we have or T~ (see Fig. [T]). For T^, C^^)'* stays +oo for Wi = +12 and Vi > —1 and — oo for Wi — —12 
and Vi < 1 whereas for T~ it is —oo for Wi — +12, Vi < —1 and +oo for Wi — —12 and Vi > 1. In summary, for any 
branch it goes all the way to — oo for negative Vi and all the way to +oo for positive Vi. Furthermore it can be made 
as negative (positive) as one wants for sufficiently large negative (positive) Vi. Thus, it is clear that independent of 
the value oi Cg^^'^ (provided it remains finite) there will always be a curve splitting the {vq^wq) plane in two parts, 
such that to the left of the curve the graviton zero mode is a ghost and to the right of the curve it is not. The case of 
the radion is more complicated so rather than giving analytic expressions describing the different behaviors, we prefer 
to show, directly in the figures, the relevant cases for the different choices of branches. 

There is also an extra source of instability in these models, namely the possible existence of tachyons in the 
spectrum. A thorough study of tachyons for the case of AdS4 branes in the standard branch was done in [§] , and the 
techniques to study them apply mostly unchanged to the self-accelerating branch. Recall that in our notation, due to 
^{^{AdSi) < 0, tachyonic solutions are characterized by q > 0. It proves useful to study the behavior of the quantity 



that has the same massive solutions as our original eigenvalue equation, but we have explicitly removed the zero mode 
solution, q = Q. The behavior of this function at g = and g ^ oo is, respectively @, 

J)(g = 0) = 2-A^Cf, (42) 

and 

sign[S)(g oo)] = sign[i;oWL(cos"^ Tq + cos"^ Tl - tt)]. (43) 

Thus, if the signs of both are different, there exists an odd number (therefore at least one) of solutions with positive 
q, i.e. tachyons. If the signs of both are the same, then there is either no tachyon or an even number of them (that 
last distinction can be resolved numerically). The condition L > gives, for AdS^ branes, Tq + > which in turn 
implies 

cos^^To + cos^^Tl -TT < 0. (44) 

It is then clear that in order to have an even number of tachyons (usually zero), we need 

sign vqVl = sign Cf^ > 0, (45) 

where the last inequality is for the phenomenologically relevant region in which the graviton zero mode is not a ghost. 

In Fig. [2] we show the case of the {Tq,T'^) branch for two possible behaviors, depending on which choice of 
parameters we make for the brane at y = L. We represent three different shades in the figure, the light one is the 
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FIG. 2: Ghost and tachyon-free region for the (Tq^jT^) branch. The different shades correspond to the following regions: 
L > (light), > (dark) and the intersection of both (intermediate). The hatched area corresponds to the region where 
> 0, L > and there are not tachyons in the spectrum. The region to the right of the thick solid line corresponds to 
' > 0. Thus, the hatched area to the right of that hne is the ghost and tachyon-free region for this configuration. In the left 
panel vl ~ 2.5, wl = 7 whereas in the right panel we use vl ~ 2.5 and wl ~ —13. 

region of parameter space for which L > 0, the dark one is the region in which the radion is not a ghost and the 
intermediate shade is the region for which both conditions are satisfied. That latter area is hatched if no tachyons 
are present in the spectrum. In these examples we have vl = 2.5 > and therefore the condition for no tachyons in 
the spectrum demands vq > 0. Finally, the thick solid line represents the points for which cf ^ = 0. Therefore the 
hatched region to the right of the solid line is ghost and tachyon free. (It can be checked numerically that the massive 
gravitons have positive kinetic coefficients in that region.) Note that there are also regions of parameter space with 
positive and small vq that are not allowed. 

In Fig. [21 we consider the case of the {T^ , Tj^) branch. In this case there is a large region where the radion is not a 
ghost for both values of the sign of uq, provided wl > 12 (considering values of wl < —12 yields a very small region 
for which the radion is not a ghost). In the figure, we have chosen vl = —4 and wl = 15. The shades have the same 
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FIG. 3: Ghost-free region for the {Tq,Tj^ ) branch, with wl > 12. We have taken vl — —4 and wl = 15. The region to the 
right of the thick solid line corresponds to Cg"' > 0. The shades correspond to the following regions: L > (light), C^, > 
(dark) and the intersection of both (intermediate). The ghost-free region has tachyons in the spectrum. 

meaning as in the previous figure. From Fig. [1] it follows that wl > 12 implies vl < for T^, thus the region to the 
right of the graviton zero mode line with vq > has an odd number of tachyons. A numerical analysis shows that the 
small region to the right of the graviton zero mode line but with vq < has two tachyons. Thus, this solution does 
not have any region that is ghost and tachyon free (no hatched area) . 

Let us consider now the cases of self-accelerating branches: {Tq ,T£) and {Tq,Tj^). In Fig. |4]we show the case 
{Tq,Tj^). There is a large ghost-free region in the case that wl > 12 whereas it shrinks to almost nothing if 
Wl < —12. From Fig. [l] it follows that wl > 12 implies vl < (left panel in the Fig. |4]) and wl < —12 implies 
Vl > (right panel). The ghost-free areas have, respectively, ug > and Uq < and therefore at least one tachyon. 
Thus in this branch there are no ghost and tachyon-free regions. In the {T^ , T^) branch, there is a large ghost-free 
region with wq > 12 (and therefore vq < —1), quite independent of the values of vl and wl- This ghost-free region. 
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FIG. 4: Ghost-free region for both the (Tq ) branch. The shades correspond to the following regions: L > (light), > 



(dark) and the intersection of both (intermediate) . The region to the right of the thick black line corresponds to Cg"' > 0. In 



the left panel we have vl 



19 whereas for the right panel «l = 4 and wl = —15. 



however, always has tachyons and therefore there are no ghost and tachyon-free regions in this case, either. We do 
not present a figure for this case since it is quite similar to the left upper region of the left panel in Fig. [J] 

In summary, in the AdSi case, only the {Tq ,T^) branch with > a-nd > Q has ghost and tachyon-free 
solutions. 

B. positive curvature: dS4, branes 

We now consider the phenomenologically more relevant case of inflating branes. This corresponds to a choice of 
brane parameters such that |Tj| > 1. Recall that the absence of a true singularity in the bulk forces us to choose 
opposite signs for Tq and T^. This leaves eight different regions of parameter space: four choices of combinations 
T^, times two choices of the sign of Tq. We have emphasized above the possible existence of a new source of 
ghosts for inflating branes, the longitudinal component of a massive graviton with mass < < 2H^. There is a 
deep connection, already present in the DGP model, between the radion and the massive graviton ghosts, that can 
be characterized by the following identity. 



S)(9 = 2)=sign(To)-(l 



vqTo){1 + vlTl)C^ 



(46) 



where we defined D in (gl]) and = kC^p/x 'H^ has the same sign as C0 for dS^. This identity states that at the 
boundary between the regions where the radion is or is not a kinetic ghost, there is always a massive graviton with 
m? = 2H^, thus also at the boundary of having a ghost longitudinal component. In the DGP model, the dependence 
on the brane parameters is such that, as we cross from the region where the radion is a ghost to the region where it is 
not, a massive graviton goes from mass squared larger than 2H^ to mass squared smaller than that, so that there is 
always a ghost in the spectrum. In our case, the richer parameter space allows for different behaviors. Let us denote 
by s = {v* ,T*) a point for which C^{s) = T){q = 2) ^ 0. The variation of the solution of the equation 



S)(g) = 0, 

around s for fixed v leads to (see Appendix [Bl for details) 



Sq_ 

6To 



8 1 dC^ 

3 B{vo,wo,vl,wl) OTo 



(47) 



(48) 



where B(vo,wo,vl,Wl) is a function that is negative in the region of parameter space for which L > and positive 
otherwise. Thus, if s is in the physical region with L > 0, the slopes of the solution of the massive graviton and the 
kinetic coefficient of the radion have opposite signs, which means that near the boundary defined by — 0, either 
the radion or a massive graviton is always a ghost. If, on the other hand the point s lies outside the physical region 
{L < 0) there is potentially a dramatic effect on the physical spectrum. In this case, contrary to what happens in 
the DGP model, there will be regions where neither the radion, nor a massive graviton is a ghost. We can study the 
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presence of massive gravitons with 5 < 2 in a way similar to the study of tachyons in the AdS^ case. The behavior of 
in the hmit q —00 reads (see Appendix |B|) 

sign[S)(g -00)] = sign[TovoVL]. (49) 

Thus, the sign of the determinant equation at minus infinity changes when either vq or goes through zero. Using 
PS)) and recaUing the definition of we also have 

sign[D(<7 = 2)] = sign[ro]eoeL, (for > 0), (50) 

where eo,L are the sign choices (branches) for T^"j^ ■ A comparison of the sign of the determinant at both points shows 
that, excluding the graviton zero mode, there is an even (possibly zero) number of modes with q < 2 {i.e. either 
ghosts or tachyons) if sign[woWL] = eo^L and an odd number of them (thus at least one) otherwise. In order to classify 
the different behaviors, it helps recalling that T[' has the opposite sign to Vi (see Fig.[l] right panel), whereas has 
the same sign as Vi everywhere except in a small wedge, where it has the opposite sign (Fig.[Tl left panel). 

Let us start our discussion with the cases for which the line = is always outside the physical region. As 
explained in Appendix [Bl this occurs in the {T^ ,T^) branch for positive To and the {T^,T£) branch with Tq 
negative or with Tq positive and (T^ — 1)(2 — vl + vlTl)/{^ + vlTl) > —4. In the latter case the radion is always a 
ghost whereas in the former two it is never a ghost in the physical region of parameter space. These are potentially 
very interesting cases, since, contrary to what happens in the DGP model, they include a self-accelerating solution 
with regions where neither the radion nor the longitudinal component of a massive graviton is a ghost. In this case 
we have sqCl = — 1 and therefore the regions that satisfy > with vqVl > have an odd number of modes with 
q < 2. This leaves the wedge of the brane in the (+) branch as the only possible ghost-free region. The massive 
graviton/radion system is actually ghost free in these regions. Unfortunately, having one of the two Vi negative makes 
the graviton zero mode a ghost. An example of this case is shown in Fig. [5l where we have chosen the {T^,T^) 
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FIG. 5: Ghost and tachyon-free region for the (Tg ,T^) branch. We have chosen vl = 0.05 and wl ~ —13. The shades 
correspond to the following regions: > (dark) and the intersection of > with L > and yo such that there are no 
singularities in the background (light). The region to the right of the thick black line corresponds to Cg"' > 0. The hatched 
region corresponds to the physical region where the radion is not a ghost and there are no modes with q < 2. In this case the 
only hatched region is to the left of the thick solid line and therefore there is no ghost and tachyon free region for this branch. 

branch with vl — 0.05 and wl = —13 (thus the brane in the (+) branch living in the small wedge). The shaded region 
of the figure represents the area where the radion is not a ghost, outside the physical region (dark shade) and inside 
the physcal region (light shade), i.e. C^, > 0, L > and yo outside the interval [0,L]. The curved lower boundary of 
the dark region represents the line Cjf, = 0, which is outside the physical region and therefore, above it, neither the 
radion nor the massive graviton is a ghost. The hatched area is the physical region in which the radion is not a ghost 
and all massive gravitons have q > 2. This region is to the left of the thick solid line and therefore the graviton zero 
mode is a ghost. 

Let us now consider the case for which the line — can be inside the physical region so that, in the neighborhood 
of such line, either the radion or the longitudinal component of a massive graviton is a ghost. This happens for the 
{Tf^,T£), {T+,Tl) branches and for the {T^ ,T+) branch for negative To or the {T+,T£) branch for positive To, 
provided {Tl - 1)(2 - vl + vlTl)/{1 + vlTl) < -4. The condition Tq • Tl < immediately tells us that there is an 
odd number of modes with q < 2 for the (Tq", T£) branch, since vq ■ < while eo^L = +1- Similarly, the (TJ*", T^) 
branch has an odd number of modes with q < 2 unless one (and only one) of the two branes lives in such small wedges. 
Numerical analysis then shows that the region with both vo,vl > is ghost and tachyon-free whereas there are two 
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FIG. 6: Ghost and tachyon-free region for the {T^,T^) branch. We have chosen vl = 0.05 and wl = —12.5. The right panel 
is just a zoomed in view of the same plot. The shades correspond to the following regions: L > and j/o such that there are no 
singularities in the background (light), > (dark) and the intersection of both (intermediate). The region to the right of 
the thick black line corresponds to Cg"' > 0. The hatched region corresponds to the physical region where the radion is not a 
ghost and there are no modes with g < 2. Thus the hatched region to the right of the thick solid line is the ghost and tachyon 
free region. 

modes with g < 2 if they are both negative. An example of this is shown in Fig. [51 where we have chosen vl = 0.05 
and Wl = —12.5 (thus the L brane in the small wedge). The color coding is the same as in the previous figure. This 
time, there is a hatched area to the right of the thick solid line (with vo > 0) that is ghost and tachyon-free. In 
the right panel we show a closer view of the region around vq = that shows an area with intermediate shade to 
the right of the thick solid line for which neither the radion nor the graviton zero mode are ghosts but there is one 
mode with q < 2 (intermediate shade unhatched region with vq < 0). For the (Tq ,TJ^) and (Tg~,TjJ') branches we 
have eoCi ~ thus again only the wedge of the brane in the (+) branch might give ghost free solutions. However, 
numerical analysis always shows two modes with g < 2 in the region where the radion is not a ghost. 

Our discussion so far has assumed a finite L. The different properties we have described naturally generalize to 
an infinite extra dimension by taking the limit L — > oo. Note however that, in the dS^ case, if j/o > Oj there is a 
different infinte extra dimension limit we can take, namely L — > yo- This corresponds to an infinite extra dimension 
in conformal coordinates, for which 

ds^ = d'^{x){-i^^dx^'dx'' + dx'^), (51) 

as shown in Appendix 1X1 The behavior is therefore different depending on the sign of ?/o- If it is negative, {i.e. if 
To < 0), then we can take the L oo limit by making = 1. This in turn is obtained for wl = 12, with vl > — 1 
if we are in the (T^) branch and vl < —1 in the {T£) branch. Now (|33l) shows in this case that the graviton zero 
mode is non-normalizable and decouples from the spectrum whereas the radion and massive gravitons remain in the 
spectrum. If j/q is positive, the limit L — > j/o is obtained by taking wl — > — oo, that makes — > — oo or vl 0, 
making Tj^ ^ — oo. In that case p4p shows that it is the radion which is non-normalizable while the graviton zero 
mode and the massive modes remain in the spectrum. Apart from the fact that either the graviton zero mode or the 
radion decouple, most of the features of the spectrum generalize to the infinite extra dimension case. One interesting 
example is the (Tq,T£) branch with negative Tq. In that case we had regions where neither the radion nor any of 
the massive gravitons was a ghost or a tachyon but the graviton zero mode was always a ghost in such regions. Since 
To < 0, however, we can take the L —^ oo limit and the graviton zero mode decouples. Thus, this is another example 
that is tachyon and ghost-free. Note that this is a standard solution but contrary to the previously found ones, it has 
VL < 0. 

To sum up, we have found that only the {T^^ ,T^) branch with both vq and vl positive is free of ghosts and tachyons 
(Fig. [6l hatched area) in the case of a finite extra dimension. For an infinite extra dimension, there are also tachyon 
and ghost-free regions with vl < in the {T^,T£) branch. We have seen examples of self-accelerating solutions, 
{Tq ,T^) with To positive for instance, where neither the radion nor any massive graviton is a ghost but the massless 
graviton is (Fig. [5l hatched area). We have seen examples of standard solutions, {T^ ,T^) for instance, where there 
is a ghost in the massive graviton/radion system (Fig. [6l light unhatched area). From these examples we conclude 
that the DGP ghost (the ghost present in the self-accelerating branch of the DGP model) is not necessarily correlated 
with the self-accelerating branch in our more general warped backgrounds. 
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V. CONCLUSIONS 

The accelerating expansion of the Universe is one of the most exciting mysteries in modern science. The existence 
of self-accelerating solutions in braneworld gravity models suggests the possibility that an extra dimension, rather 
than dark energy, may explain this phenomenon. Unfortunately the original DGP model has been criticized for the 
presence of ghosts in the weakly-coupled regime of the self-accelerating solution. Generically either the radion or the 
longitudinal component of a massive graviton is a ghost. 

Given the obvious importance of self-accelerating solutions, we have investigated the presence of ghosts in a gen- 
eralization of the DGP model. Our model uses a warped AdS^ bulk and two branes with arbitrary tensions and 
brane-localized curvature terms. Depending on functions of the input brane parameters which we have denoted as 
To, Tl, the AdS5 bulk is sliced into sections of negative, positive, or zero curvature, giving rise to AdSi, dSi, or M4 
branes, respectively. For each value of the input brane parameters, there are in general two different static back- 
grounds per brane; the four resulting branches are denoted by {T^,T^). In the dS^ or flat case, the (Tq~,T^) and 
(Tf^,T£) branches provide self-accelerating solutions. 

We have performed a comprehensive analysis of the spectrum for both AdS^ and dS4 branes and every branch, with 
special emphasis on the presence of ghosts. The results can be summarized as follows. For the ^^54 brane on the 
{T^,T^) branch, the spectrum is free of ghosts and tachyons provided both localized curvature terms are positive; 
on all other AdS4 branches, for all values of the input parameters, the spectrum has either ghosts or tachyons. For 
^5*4 branes on the {T^ , T^) branch the spectrum is also free of ghosts and tachyons provided both localized curvature 
terms are positive. On all other dS4 branches, with finite L, for all values of the input parameters, the spectrum has 
either ghosts or tachyons. This feature of our results is discouraging. 

We find, however, solutions that depart from the general behavior of the DGP model in interesting ways. We have 
self-accelerating solutions which are free of DGP ghosts, i.e., neither the radion nor any massive graviton are ghosts; 
in these cases the graviton zero mode turned out to be a ghost. In some of these solutions, one can take the second 
brane to infinity so that the (ghost) graviton zero mode decouples, thus obtaining another example of ghost and 
tachyon-free model. We also have solutions which are not self-accelerating, but nevertheless have a DGP ghost. 

Although there is no self-accelerating solution that is totally ghost-free, the fact that the ghost present in the 
DGP model is not correlated with the self-accelerating branch makes these models an excellent laboratory to study 
self-acceleration. One can hope that simple modifications of these models could get rid of ghosts altogether. For 
example, in the case of a flat bulk it was recently suggested that adding a Gauss-Bonnet term to the bulk action 
will exorcise the DGP ghost [2^. Since in some cases the ghost is a radion, one could also attempt to modify the 
bulk scalar dynamics, in the style of Goldberger and Wise '27] (see also [23]). In this context, it is important to 
also consider the effects of brane matter, which requires a fully time-dependent analysis rather than the quasi-static 
approach taken here. Another promising direction is the study of exact solutions, using the methods developed in 
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APPENDIX A: COORDINATE SYSTEMS 

In this Appendix we will give the AdS5/{A)dS4 background solutions in a different coordinate system that has 
been used in other studies of brane induced gravity models and is also the preferred one in cosmological studies. This 
coordinate system is conformal: 

ds^ = d^{x){'yf,^dx''dx'' + dx^), (Al) 

where 

a^ly) = cL^lx), dy^ = d'^{x)dx^. (A2) 
The coordinate change is of course different for the dS4 and AdS4 cases. 
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• AdS4 case: 

n ( = 

k smH{x — xq) 



«(^) = T-:x-E77^-rT' (A3) 



with 



1 



/sinhA;(y — t/o) — i 



and 



y = yo-T sinh ^ [ cot H{x - a;o)] . (A5) 



Requiring that a; = for ?/ = we obtain 

2H \ sinh kyo + i J 



i /sinhfcyo-i\ 



In particular it is easy to check that 



y-yo^O^ => x-xo^ttt/{2H), (A7) 
?/— >oo => a; — ajQ— >0~. (A8) 



• dS4 case: 

ninfi — 

k sinh if (a; — a;o) ' 

with 



aW = f „:„.J„ (A9) 



1 /coshA;(y-yo) - 1\ 
"-^° = -2g^H coshfc(.-,o) + l J' 

and 

y = yo — ^ cosh~^ [cothH{x — xo)], (All) 

where eo = sign Tq. Again, requiring that x = for y = we obtain 

1 , / cosh kyn — 1 \ / . . „^ 

Xo= n — . (A12) 

In this case we have 

y^yo =^ X ^ oo, (A13) 

y-*oo ^ X- xo ^ 0- (A14) 

APPENDIX B: DETAILED CALCULATIONS FOR dS4 BRANES 

We reproduce in this appendix the detailed calculations that lead to the determination of ghosts and tachyons for 
dS4 branes. For convenience, we will use {v, T) as independent variables. The first step is to characterize the region 
of parameter space satisfying (7^ = 0. Let us fix the value of {vl,Tl), and solve C^{vq,Tq) = for VqI 

-2T,*-2-XL 



where 



YT^L • ^^'^ 
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In the (uq, '7o)-plane, the allowed regions are bounded by Tq = ±1 or Tq = —Tl, and Tq = —I/vq. For example, the 
allowed region for the (Tq,T£) branch with < lies in the fourth quadrant with Tq < —T^ and Tq < —1/vq. In 
this case, in order for the = line to fit in the allowed region, it should be satisfied that Vq of (jBl|l be larger than 
— 1/Tg*. This inequality can change sign for some of the branches. 

Therefore, to see whether (7,/, = line lies inside or outside of the allowed region, we need to check the sign of 



whose only nontrivial part is 



1 

n 



1 



1 



T* T*2 + (2 + xl)T^ + 1 



(B3) 



e(To) =T*^ + (2 + xl)To* + 1. 



(B4) 



Since Tq < — Tl < — I or 1 < Tq < ^Tl, from the elementary analysis on quadratic functions it would suffice to 
determine the signs of 



6(-Tr 



1 



1 + vlTl 



6(1) ^4 + XL, 



(B5) 



and the location of the symmetry axis, I — —(2 + Xl)/2, relative to 
opposite of €l- The result is summarized in the following table: 



-Tl or 1. Note that the sign of &{-Tl) is the 





&{-Tl) e(i) I 


sign[(lB3|] 


allowed region 


location of = 


To < 






H — > — as To increases 


VQ + 1/To < 


inside or outside 




+ > -Tl 


+ 


VQ + 1/To < 


outside 


iTo,Tl) 




H !■ — as To increases 


VQ + 1/To > 


inside or outside 


iTo,TL) 


+ > -Tl 


+ 


VQ + 1/To > 


inside 


To>0 






+ 


VQ + 1/To > 


inside 




+ + < 1 
or — 


H > — as To increases 


VQ + 1/To > 
VQ + 1/To > 


outside 
inside or outside 






+ 


VQ + 1/To < 


outside 


{To,TZ) 


+ + < 1 
or — 


^ — > — as To increases 


VQ + 1/To < 
VQ + 1/To < 


inside 
inside or outside 



A blank in the tables denotes that the corresponding information does not affect the result. H > — as To increases 

means that both signs are possible, depending on the parameters. Comparing the sign of (|B3[) with the allowed region 
for each case (columns 4 and 5 in the table) we obtain the classification of the last column in the table. 

Let us assume the line of = is inside the allowed region, with {vq,Tq) representing the locus. The treatment 
of the solution is different depending on the sign of Tq. If Tq is negative, so that z is positive, we get, from (|30l34p . 



S)((7 = 2) = 4(To + 1)(2 + VQ + vqTq){1 + vlTl) + 4(Tl - 1)(2 - vl + vlTl){1 + vqTq) 
= -^{l + vqTq){1 + vlTl)C^ , 



(B6) 



and therefore if = has a solution, there is a KK graviton mode with q = 2. Let us now consider how and 
J) change as we vary To around the solution: for given L-parameters, when s = {q = 2,Tq = Tq) solves J) = 0, at 
(q — 2 + Sq, Tq — Tq + STq) which is another solution in a small neighborhood of s, we have 



= D{T* + 6TQ,2 + Sq) = D{s) 



STq 



= 0-^il + VoTQ){l + VLTL)^^ 



93 

dq 



Sq 



STq + — 



5q, 



(B7) 



I.e. 



5q_ 

STq 



{1+VqTq){1 



d^/dq\ 



VlTl) dC^ 



(B8) 
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To evaluate dD/dq, we use ^ 3.7 (6), fsq 4.224.9 and [3i| 3.6.1 (8) to get 



5.P°(z)|.=o 
a.P-^(z)|.=o 



In- 



z- 1 
z + 1 
3^ + 1 



z + 1 , z + 1 
In ■ 



z - 1 

2 + 1 



4(z + l) 2(z-l) 
Also using [30*1 8.712, 2.727.2 and [3l| 3.6.1 (5) we obtain 

5.Q°(z)|.=o 



■In- 



2 ' 

z + \ 



iln21n^ + iln(z2_l)ln^ 
2 z- 1 4 ^ ' z+1 



9.Qi(z)|.=o 



2^/ 



- 1 



; In ■ 



z + 1 
z - 1 



In- 



— z 



2(z2- 1) 

where Li2(z) is the dilogarithm function, defined as 



1 f / 9 1 M -2 + 1 



2zln 



z2- 1 



fc=i "^^ 



° ln(l - t) 



dt. 



Then, dT)/dq\q^2 becomes 



dq 



4(1 + VLTL)[2{n + 1)(1 + i'oTo) + (2 - 1-0 + voTo){-n + 1) In ^ } 



9=2 



+ (2 - + vlTl){Tl - l){(«o + 2To + v^T^) In ^ | + 2(1 + v^Tq) (l + In ^^-^) } 



T2-1- 



+4(1 + vaTo)[2{TL - 1)(1 + i'lTl) - (2 + + i;LTL)(rL + 1) In ^^^^} 



^(2 + wo + t'oTo)(To + 1){{VL + 2ri + vlTD In - 2(1 + vlTl) (l + In ^^^) } 



In - ^ 



At To = To* , (|B6l)= simplifies the above into 



as 

9g 



(To*2 - l){Tl - 1) 8(r* + Tl) - xl In + (xl + 4)^ In 



and (IB8I) becomes 



where 



(Tl - 1)2 - xlTl 



Sq_ 

6To 



[T* + 1)2 + xlTJ 



3B(To; Ti,XL) dTo 



B = S{T* + Tl)-xI In- 



'Tn* - 1 



" ^+(XL+4)2ln- 



-B is a growing function of Tq, 



Tl-1 



dB _ (2 + 27o + xl)^ 
dTo"' T2-1 



Tl + 1 



(B9) 
(BIO) 
(Bll) 



(B12) 
(B13) 
(B14) 

(B15) 



(B16) 



(B17) 



(B18) 



(B19) 



(B20) 



For L > we have Tq < —Tl and therefore B is bounded above by B{—Tl) = 0, or in other words, it is always 
negative. If L < 0, on the other hand, we have Tq > —Tl and therefore B is bounded below by 0, so it is positive. 
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Thus, when = has a solution in the physical region (with L > 0), the slope of the mass of the first KK graviton 
mode when it passes g = 2 and that of when it crosses zero have opposite signs, and therefore when the first KK 
graviton mode is heavier than 2H^ the radion is a ghost and vice versa. If = has a solution outside the physical 
region (where L < 0), then both the radion and the longitudinal component of a massive graviton are ghosts at one 
side of the line (in the {vq,wq) plane) and neither of them is at the other side. 

If To is positive, then z is negative in the physical region and we have to use the wave functions evaluated at —z 
for the massive gravitons. This introduces some extra minus signs, so that now we have 

S(g = 2) = ^{1 + voTo)il + vlTl)C^ , (B21) 



Sq_ ^ 8 {1 + vaTo){l + vlTl) dC^ 

STo s 3 dD/dq\s dTo 

as ^ (r*^-i)(r|-i) sjr^ + n) - xl in + (xl + 4)^ in 

dq s in - 1)2 - xlTl ■ (T* + 1)2 + xlTo* 



(B22) 
(B23) 



These extra minuses get cancelled among themselves, and we still get the same final result, (|B18IIB19l) . Thus inde- 
pendently of the sign of Tg, it is always true that, near the locus of points satisfying = 0, either the radion or 
the longitudinal component of a massive graviton with < 2H^ is a ghost if parameters are chosen from inside the 
physical region and neither of them is to one side of the region = if that region has L < Q. 

Let us now compute the limit D{q — oo), which is equivalent to v ^ oo, where is the index of the corresponding 
Legendre function. The calculation is again different for the two signs of Tq. We start considering the case of negative 
To, so that z is positive. Using 8.723 1 from [33| we find 

1 p(i^+l)a 

p-i^(cosh«)~^^===i.-5/2, (B24) 
where we have defined z = cosh a. Then 

(1 - z')^p-^^{cosha) ^ --}-^^f^^Xe^-y-^l\ (B25) 
Using now 8.723 2 and 8.732 1 from [30(, respectively, we get 

p-va 

g2_^(cosha) ~ V^_==z.3/2^ (B26) 
Ve^° — 1 

(l-^^)^Qf,^,o(cosha) ^ ^y^5^e^(^+i)"^.5/2. (B27) 

These expressions allow us to compute the limit of the different coefficients, that turns out to be identical for both 
branes 

a.(g--oo) ^ "'^ (B28) 
hM - -oo) ^ ^v,{Tl - 1) /^'J' (B29) 
We can now compute the limit of the full determinant, 



^^0-1 7^1-1 3 

^(P-ota _ 1 ^e2"i - 1 



Ve2"" - 1 Ve2"^ - 1 ^ ^ 

where we have used that, in the physical region with negative Tq, we have < T^, < — Tq ^ cosha^ < coshao 
OLL < ctQ- Thus, 

sign[S)(5 -00)] = -sign[i;ofL] = sign[roWoWL] ■ (B31) 
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We can now check what happens in the case that Tg > so that z < 0. As we said above, the vaUd solution is then 
the corresponding associated Legendre functions evaluated at — z. The expression of the and bi is then identical to 
(f30| with the change 9i —9i, i.e. 



JTo>0) 
^0 



(To) 



(B32) 



where Ci stands for any of a.; or bi. Thus, noting that now < Tq < — 7l coshap < cosha^ =J> < Q^l in the 
physical region, we get 



Tn>0 



' Ve^^o - 1 Ve^^i- - 1 

n - 1 r2 - 1 



-VqVl- 



and therefore 



' Ve^^o - 1 Ve^"^ - 1 
sign[£>(^«>''^((7 -oo)] = sign[wo?;L] = sign[ToUoWL] ■ 



(B33) 



(B34) 



APPENDIX C: FLAT BULK LIMIT 

In this appendix we will show how to obtain the limit of a flat bulk (fc 0) for dS'4 branes. Taking the fc — > 
limit in the AdS5/dS4 background needs special care because our usual reparametrization of the extra dimensional 
coordinate y, z — cothfc(y — i/q), breaks down and the input parameters, Vi = kMf /M^ and Wi = Vi/2kM^, are ill- 
defined. Therefore, it is easier to redo our analysis directly with k = 0. That is, we work on an M5/dS4 background. 
Let us define the following combinations of brane parameters and the fundamental Planck mass, 

A-^^ U - (CD 
Solving for the background is straightforward, and we get 

Gfj^j^dx^'dx^ = a2(y)7^,dx^da:'' + dy^ , (C2) 

where a = 1 + eoHy with eg = +1 for the self-accelerating branch and eo = — 1 for the standard one. Here 7^^, is the 
dS4 metric. The values of H and L are determined by the brane-boundary equations: 

= ^K=-.„H-i|, (C3) 

We first look at the radion. Following footsteps of §5 of [§], we can get the equations of motion for the massless 
scalar degrees of freedom. Solving them, we obtain 

/ N f2ix) C(x) D(x) ^, , ,, , 

V2{x,y) = h{x) + [^T[y)+'H^:S{y))^{x). (C6) 
f2{x) turns out to be a pure gauge field and will be killed. Plugging (|C5[) into the brane-boundary equations gives 

a,C + P^D ^ , (C7) 

At 



where 



i(3eo - l)/2 • a- + Xitia^ ^^^^ 



a,; 



^ g.(36o + l)/2 ■ a-(^+-o)/^ + A.Ha-^^+-")/^ ^ ^^^^ 
9i + eaXiTi/ ai 



20 



with Oi = a(yj). 

At the massive level, a solution of the bulk equations of motion is now 



b,^{x, y) = a(i-^)/2A,.(x) + a(i+^)/2i3^.(a;) , (CIO) 
and the brane-boundary equations provide the determinant equation, which gives the mass spectrum; 

= {(3+ V9^)eo-WAo(?}-a(f'+^/^)/'{(-3+ v/9^)(eo+WL)-WAig} 

-(3+v'9^)/2 



- {(-3 + ^9^)60 + nXoq} ■ al^'+^'-^''>^'{{S + ^/9~M)ieo + HL) + TiA^g} . (Cll) 

Let's look into the self-accelerating branch(eo = +1). Here, since ao = q;l = ^1, we can fix the gauge function 
F{y) such that J'{yi) = xPi with the constant x ^ residual gauge parameter. This choice kills the C mode and 



D^^^, (C12) 



which redefines (pS)) for M5 / dS4 background; 



yi{y) = ^a~' - d , y2{y) = xUa-^ + . (C13) 
ri a 

Working out the quadratic action, we finally get 

= \- E {y^ + \/H ~ 9. + H{y, + A.)) ■ ^^^^^ 

Also with eo — +1, (|C11|) becomes 

= {(3+ v/9^)-7iAo<7}-air^+^^^^/'{(-3+y9^)(l + WL)-7iALg} 

- {(-3 + ^9^) + nX^q} ■ a^<^+^/^'/'{(3 + v/9^)(l + HL) + HAl?} . (C15) 

We can match the above with the recent results of Izumi et al. for braneworld models on a M5 /dS4 background with 
a finite extra dimension 28]. The conversion from our setup to theirs can be achieved by the following identifications: 
X{),Xl 2rc, Vi Ti. Then correspondence among derived quantities follows: Ti = + l/rc L 

— q mf. Upon making these substitutions, (|C14[) and (jClSP turn into 

= 3x^(^,-i/_^l + (l + ^^+^^)(^ + ^^-^^) 



(l + 2i/+rc)(l + 2H_rc) 

This expression is a positive constant times times the expression in eqn. (2.42) of [2^. Similarly: 



= {3-2H+ramf + Jg-Am^ 



i7+\{-l+V9-4™?)/2 



3 - 2H_rcmf + \/9- Am) 



- (3 - 2H+rcmf - ^9 - 4TOf ) 



^ ^ (-3-2i/_r,mf-v'9-4mf), (C17) 



X 



which is 4H^ times the determinant of the expression in eqn. (2.35) of f2§\. The normal branch solutions (ep = — 1) 
can be worked out similarly. 
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